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A  COMPREHENSIVE  MATHEMATICAL  MODEL  OF  THE  CARDIOVASCULAR  SYSTEM 
UNDER  TIME-DEPENDENT  ACCELERATION  STRESS 


Xavier  J.  R.  Avula 
Department  of  Engineering  Mechanics 
University  of  Missouri -Roll a 
Rolla,  Missouri  65401 


ABSTRACT 


In  this  study  a  comprehensive  mathematical  model  of  the 
cardiovascular  system  under  time-dependent  accelerations  is  developed. 
Recently  developed  high  performance  aircraft  would  expose  the  human 
body  to  acceleration  injury  if  appropriate  life-supporting  devices  are 
not  incorporated  in  the  design.  To  aid  in  the  construction  of  desirable 
life  support  systems  for  aerospace  maneuvers,  the  deformation  of  the 
arterial  and  venous  segments  under  dynamic  fluid  loads  caused  by  blood 
pooling  during  Gz  acceleration  are  calculated.  Linearized  Navier-Stokes 
equations  for  blood  flow  and  equations  of  large  elastic  deformation  theory 
for  blood  vessel  deformations  are  used.  The  resulting  nonlinear  partial 
differential  equations  are  solved  numerically.  The  model  presented  here 
consists  of  a  closed-loop  hydrodynamic  system  including  the  heart  pump, 
compartments  of  large  arteries  and  veins  in  the  upper  and  lower  body,  and 
a  baroreceptor  feed  back  mechanism.  To  verify  the  model  aortic  pressure 
is  calculated  for  an  experimental  deceleration  profile.  A  satisfactory 
agreement  between  the  theory  and  experiment  is  found. 

NOMENCLATURE 


A  =  Material  constant  in  the  strain  energy  function 


B  *  Material  constant  in  the  strain  energy  function 


B  *  (with  sub-  and  superscripts)  *  Surface  metric  tensor 
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Information  Division 


F  *  Body  force 
f  *  Acceleration 

6  =  Determinant  of  the  covariant  metric  tensor  6^ 
g  =  Determinant  of  the  covariant  metric  tensor  g. . 

*  J 

G. .  *  Covariant  metric  tensor  for  reference  coordinates  of  the 
1J  deformed  state 

G1J  *  Contravarient  metric  tens or  for  reference  coordinates  of 
the  deformed  state 

g. .  =  Covariant  metric  tensor  for  reference  coordinates  of  the 
undeformed  state 

g1J  =  Contravarlant  metric  tensor  for  reference  coordinates  of 
the  undeformed  state 

2h  =  Wall  thickness  of  the  blood  vessel 
o 

*1**2’*3  =  Strain  Invariants 

L  =  Length  of  blood  vessel  element 
n°^  =  Stress  resultant  tensor 
p  =  pressure 
Q  *  Flow  rate 

=  Internal  radius  of  the  blood  vessel 

r  =  Radial  coordinate 
t  =  time 

u,w  =  Radial  and  axial  fluid  velocity  components 
W  =  Strain  energy  function 
z  =  Axial  coordinate 

r“  =  Christoffel  symbol  of  the  second  kind 
pp 

Y  =  Shear  rate  In  blood 
A  =  Stretch  ratio 

p  =  Coefficient  of  viscosity  of  blood 
v  «  Kinematic  viscosity  of  blood 
p  =  Mass  density  of  blood  vessel 
pQ  *  Mass  density  of  blood 

t1'-’  =  Stress  tensor 

4>  =  A  function  of  strain  energy 

f  *  A  function  of  strain  energy 


Subscripts  and  Superscripts 
i,  j,  k  =  1,  2,  3 


r,  s  =  1,  2 

a,  8  -  1,  2 

INTRODUCTION 

The  human  body  is  well  accustomed  to  the  earth's  force  of 

gravity,  but  recent  space  age  developments  have  occasioned  its  exposure 

to  the  hazards  of  high  and  abnormal  gravitational  fields,  which  are 

manifested  in  the  form  of  vibration,  impact,  weightlessness,  and  positive, 

negative,  forward,  backward,  and  angular  accelerations  that  are  beyond 

its  tolerance  levels.  Abnormal  accelerations  on  the  human  body  are  known 

to  cause  headache,  abdominal  pain,  impairment  of  vision,  hemorrhage,  and 

1-3 

fracture  depending  upon  the  severity  and  kind  .  The  cardiovascular 
system  in  being  central  to  the  homeostasis  of  the  organism  is  extremely 
susceptible  to  hostile  changes  in  environmental  force.  The  design  of 
protective  devices,  which  are  expected  to  provide  acceleration  tolerance 
for  the  organism  during  aircraft  and  spacecraft  maneuvers,  must  take  into 
consideration  the  response  of  the  cardiovascular  system  to  acceleration 
stress.  Therefore,  a  thorough  understanding  of  the  system  and  Its  structure- 
function  relationship  in  an  abnormal  force  environment  is  essential  to  any 
effort  directed  to  overcome  the  acceleration  trauma. 

The  prohibitiveness  of  actually  subjecting  the  human  body  to  abnormal 
accelerations  to  gain  knowledge  of  the  cardiovascular  system's  response  is 
obvious.  The  alternative  Is  to  develop  a  mathematical  model  and  to  investigate 
the  response  of  the  system.  The  need  for  mathematical  models  and  the  analysis 
of  model  features  for  prediction  of  system  performance  are  well  recognized 
In  view  of  the  cost  and  risk  Involved  In  testing  the  original  system. 
Theoretical  analyses  are  extremely  helpful  for  evaluating  the  relative 


injury  potential  for  various  acceleration  functions,  in  guiding  experi¬ 
mental  investigations,  and  in  developing  and  understanding  protective 
measures.  Mathematical  procedures  also  provide  the  basis  for  establishing 
precise  dynamic  and  physiological  scaling  laws  needed  to  translate  experi¬ 
mental  data  obtained  with  various  species  into  meaningful  results  for 
humans . 

There  is  no  dearth  of  mathematical  models  of  the  cardiovascular 

4  5 

system  in  the  scientific  literature.  Womersley  and  Noordergraaf  presented 

a  mathematical  analysis  of  some  aspects  of  the  cardiovascular  system  by  using 

6  7  8 

a  lumped  parameter  model.  Taylor,  Kenner,  and  Attinger  et  al.  used 

distributive  parameter  models  to  analyze  pressure  -  flow  relationships  in 

arteries  and  veins.  It  is  generally  believed  that  a  lumped  parameter  model 

is  superior  to  the  distributive  one  in  the  evaluation  of  overall  cardio- 

g 

vascular  system  performance.  Beneken  and  DeWit  characterized  a  large 
analog  model  of  the  entire  human  circulatory  system  in  the  form  of  approxi¬ 
mately  40  equations.  Rigorous  mathematical  formulations  and  extensive 
physiological  information,  including  the  factors  affecting  contractil ity 
of  the  myocardium,  the  effect  of  intrathoracic  and  abdominal  pressure 
changes  on  the  venous  conductance  and  the  ventricular  distention,  the 
baroreceptor  reflex  control  of  the  cardiac  pump  and  vascular  resistance, 
the  autoregulation  in  various  vascular  beds,  and  capillary  fluid  shift 
and  stress  relaxation  phenomena  were  introduced  into  the  model.  It  has 
been  demonstrated  that  the  model  provides  quick  solutions  for  parametric 
sensitivity  tests.  Guyton  and  Coleman1®  presented  an  analog  model  of 
long-term  circulatory  regulation  and  emphasized  the  integration  of  the 
long-term  autoregulation  of  the  systemic  vascular  system  into  the  basic 
scheme  as  the  most  powerful  control  mechanism  for  tissue  homeostasis. 

Pennock  and  Attinger11  proposed  a  mathematical  model  to  analyze  the 
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overall  performance  of  the  oxygen  transport  system.  This  model  was 

represented  by  six  equations  and  used  to  describe  oxygen  transfer  and 

transport  and  normal  performance.  The  changes  on  the  interaction  of 

subsystems  and  on  optimization  with  respect  to  flexibility  and  maximal 

limits  of  performance  were  examined  by  altering  various  parameters. 

McLeod  proposed  a  physiological  simulation  benchmark  experiment  (PHYSBE) 

to  economize  on  programming  efforts  and  to  establish  bases  for  a  comparison 

13  14 

of  various  types  of  computer  models  of  circulation.  Sagawa  *  described 

the  overall  circulatory  regulation  and  the  mechanical  properties  of  the 

15 

cardiovascular  system  in  the  control  of  circulation.  Camill  studied 

the  response  of  the  human  cardiovascular  system  to  whole-body  sinusoidal 

vibrations  by  using  an  open-loop  analog  model  of  the  standing  and  sitting 

man.  Boyers  et  al.^  simulated  the  steady  state  response  of  the  human 

cardiovascular  system  with  normal  responses  to  change  the  posture,  blood 

loss,  transfusion,  and  autonomic  blockage.  Collins  et  al.17  presented  a 

dynamic,  mathematical  simulation  of  the  cardiopulmonary  system.  Several 

articles  related  to  blood  flow  in  arteries  have  appeared  in  the  book  by 
18 

McDonald  .  An  elastic  tube  theory  of  blood  flow  has  been  treated  by 

iq  20  21 

Lambert  and  Skalak  and  Stathis  .  Kivity  and  Collins  presented  a 

viscoelastic  tube  model  for  aortic  rupture  under  decelerative  forces. 

22 

Rudinger  studied  the  effect  of  shock  waves  on  mathematical  models  of 
aorta  for  better  understanding  of  the  behavior  of  the  actual  aorta. 

Most  of  the  above  modeling  efforts  deal  with  electrical  analogs 
of  the  cardiovascular  system  In  which  various  parameters  are  introduced  in 
terms  of  resistances,  Impedances,  and  capacitances.  Because  the  measured  or 
postulated  values  of  these  parameters  are  variable  and  their  representation 
under  high-g  conditions  Is  speculative,  errors  of  large  magnitude  are 
likely  to  creep  Into  the  model.  Therefore,  an  analysis  that  is  based 
purely  on  the  original  properties  of  the  cardiovascular  elements  coupling 
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the  fluid-flow,  deformation  of  the  vascular  walls,  and  the  material 

properties  of  the  blood  vessels  and  surrounding  tissue  is  expected  to 

yield  a  better  mathematical  model. 

To  understand  the  blood  flow  characteristics  in  the  arterial 

system,  the  knowledge  of  the  material  properties  of  the  arterial  wall 

23  24  25 

is  essential.  Bergel  ,  Fung  ,  Demi  ray  and  Vito  have  utilized  mathe¬ 
matical  models  of  the  constitutive  properties  of  the  arterial  tissue  to 

determine  the  stresses  in  the  arterial  walls.  In  the  present  study,  the 

25 

strain  energy  function  given  by  Demiray  and  Vito  for  an  arterial  wall 
specimen  has  been  used  in  determining  the  aortic  pressure  that  is  compat¬ 
ible  with  large  deformation  of  the  aorta  and  the  associated  flow  under 
acceleration  stress. 

The  present  study  also  considers  the  effect  of  acceleration  on  the 

mircrocirculation.  Microcirculation  under  normal  conditions  was  investigated 

26  27  2ft 

by  Prothero  and  Burton  ,  Whitmore  ,  Gross  and  Aroesty  ,  Gross  and 

29  30  31 

Intaglietta  ,  Skalak  and  Fung  who  presented  various  theories  of  flow 

in  the  capillary  bed  connecting  the  arteries  and  veins. 

Several  experimental  investigations  on  the  effects  of  acceleration 

stress  on  the  human  body  have  been  performed  at  the  USAF  School  of  Aerospace 

32 

Medicine  at  Brooks  Air  Force  Base,  Texas.  Burton  subjected  miniature 

swine  to  Gz  acceleration  to  study  its  effects  on  the  organism  and  extra- 

33 

polated  the  results  to  human  beings.  Parkhurst,  et  al  conducted  experi- 

ments  on  human  tolerance  to  high  +GZ  forces.  Leverett,  et  al.  investigated 

the  physiologic  response  to  high  sustained  acceleration  stress.  Peterson, 

35 

et  al.  studied  the  cardiovascular  responses  during  and  following  exposure 

to  +GZ  forces  In  chronically  Instrumented  anesthetized  dogs.  Burton  and 
36 

MacKenzIe  determined  the  extent  of  heart  pathology  as  a  function  of 
acceleration  stress. 
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DESCRIPTION  OF  THE  PHYSICAL  MODEL 

Because  the  cardiovascular  system  consists  of  several  components. 

It  would  be  too  complicated  to  handle  all  of  them  In  a  general  model. 
Although  It  would  be  desirable  to  include  the  behavior  of  each  cardio¬ 
vascular  component  under  acceleration  stress  in  the  total  modeling  effort, 
certain  components  can  be  lumped  together  to  simplify  the  analysis  and 
still  preserve  the  character  of  the  system.  In  the  proposed  model,  five 
elastic  chambers  containing  blood  are  arranged  in  a  closed  loop,  and  a 
mathematical  analysis  is  made  to  calculate  the  fluid  shift  under  acceler¬ 
ation  stress.  The  five  elastic  chambers  are:  1)  arteries  of  the  thorax 
and  the  lower  body,  2)  veins  of  the  thorax  and  the  lower  body,  3)  arteries 
of  the  upper  body,  4)  veins  of  the  upper  body,  and  5)  heart  and  lungs. 

A  schematic  of  the  elastic  reservoirs  arranged  in  a  closed  loop  is  shown 
in  Fig.  1.  This  physical  system  is  then  subjected  to  +GZ  acceleration. 

The  elastic  reservoirs  are  considered  highly  deformable,  and  the  theory 
of  large  elastic  deformations  is  applied  to  the  calculation  of  their 
expansion  under  acceleration  stress.  The  Navier-Stokes  equations  are 
used  to  determine  the  fluid  velocity  into  and  out  of  the  chambers  with 
proper  boundary  conditions  to  match  the  wall  motion. 

MATHEMATICAL  MODELING 
A.  Equations  of  Fluid  Motion 

The  geometry  of  the  elastic  tube  containing  blood  in  motion  Is 
shown  In  Fig.  2.  Let  r,  e,  z  be  the  cylindrical  polar  coordinates  and 
let  u,  v,  and  w  bet  the  velocity  components  in  the  corresponding  directions. 
Assuming  axial  symmetry  In  flow  and  tube  deformation,  the  linearized  Navier- 
Stokes  equations  for  the  flow  of  blood  can  be  written  as: 


-7- 


CAPILLARY  BED 


ARTERIES  IN  THE 
UPPER  BODY 


ARTERIES  IN  THE 
THORAX  AND THE 
LOWER  BODY 


CAPILLARY  BED 


VEINS  IN  THE 
UPPER  BODY 


VEINS  IN  THE  THORAX 
AND  THE  LOWER  BODY 


Fig.  1  Schematic  of  physical  model 
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where  p  is  the  pressure,  v  is  the  kinematic  viscosity,  pQ  is  density  of 
blood  and  g(t)  is  the  body  force  per  unit  mass  caused  by  the  acceleration. 
The  continuity  equation  is 


|u  +  u+|w  =0 
3r  r  3z 


(3) 


The  above  equations  are  nondimensionalized  using  a  typical  length,  Rq, 
which  is  the  initial  (undeformed)  radius  of  the  aorta,  and  U,  the  average 
velocity  of  blood  in  the  aorta.  Introducing  the  new  quantities 


the  equations  of  motion  and  the  continuity  equation  in  terms  of  the  newly 
defined  variables  become 
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Deleting  the  "stars"  for  simplicity,  the  governing  equations  in  the 
dimensionless  form  will  become 
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The  boundary  and  initial  conditions  are 
dR, 


u  s  dT  at  r  ’  Ri 


t  >  0 


w  =  0  at  r  =  R, 


t  >  0 


(ID 


w  =  1  at  z  =  0 


t  >  0 


where  R^  is  the  inside  radius  of  the  blood  vessel  in  the  deformed  state. 

B.  Equations  of  Motion  for  Thin-Walled  Elastic  Tube: 

The  theory  of  large  elastic  deformations  is  utilized  to  describe 

the  time-dependent  deformation  of  the  blood  vessels.  In  view  of  the 

published  results  on  blood  pooling  and  the  consequent  cardiac  insufficiency, 

the  application  of  large  deformation  theory  appears  necessary.  Demiray  and 
25 

Vito  have  previously  used  this  theory  to  calculate  the  deformation  of 
arteries. 

The  undeformed  and  deformed  cylindrical  tubes  are  shown  in  Fig.  3. 

Let  r,  e,  z  represent  a  point  in  the  wall  of  the  undeformed  tube,  and  R,  0, 
z  in  the  deformed  tube,  rj,  r2  are  inside  and  outside  radii,  respectively, 
of  the  undeformed  tube,  and  Rj,  R2  those  of  the  deformed  tube.  Axial  stretch 
of  the  tube  is  neglected  because  of  tethering  caused  by  the  surrounding  tissue 
Assuming  the  material  of  the  blood  vessels  to  be  homogeneous,  incompressible, 
and  isotropic,  the  stress  at  any  point  can  be  written  as: 

t*'*  =  <f>g^  +  y/B1'*  +  PG1^  (12) 


-10- 


where  *  =  2(3W/ai1),  <p  =  2(3W/3l2),  Blj  =  I1gij  -  girgjSGrs.  Pisa 

scalar  function  which  represents  a  hydrostatic  pressure,  W  is  the  strain 

energy  function,  1^  and  1^  are  the  strain  invariants,  and  g’d,  g^,  G^, 

37  38 

and  G^.  are  the  contravariant  and  covariant  metric  tensors  *  .  The 

indices  i  and  j  take  the  values  1,  2,  and  3.  The  equations  of  motion 
are  given  by: 


T1J  I  •  +  P 
1  1  ^W 


Fd  =  p,.,f j 


(13) 


where  ||  denotes  covariant  differentiation,  pw  is  the  density  of  the 
vessel  wall,  F  is  the  body  force,  and  f  is  the  acceleration.  Let  us 
neglect  the  body  force  on  the  vessel  wall  in  comparison  to  its  effect 
on  the  fluid  flowing  in  the  cylindrical  tube.  Performing  the  covariant 
differentiation  on  the  remaining  part  of  the  equation  of  motion  we  get 


4.  r1  S3  x  r1  S*  =  n  fJ 

Hw 


T  »1  +  r  ir  T  +  r  ir  T 


(14) 


i  37  38 

where  r  ^  represent  the  Christoffel  symbols  of  the  second  kind  ’  . 

It  has  been  shown  that  for  a  biomaterial,  a  reasonable  strain  energy 
function  as  shown  in  Ref.  25  is 


o  A(I«  ■  3) 

“  ■  Ir  t*  2  -« 


(15) 


in  which  A  and  B  are  material  constants.  Defining  the  circumferential 
stretch  ratio  X  =  R/r,  the  stresses  in  the  r,  9,  z  directions  can  be 
expressed  as 

11  1  a(I2”35 

t  1  *  P  +  B(1  +  ft)  e  2 


i  ii  i  A(I?“3) 

=  P  +  B(1  +  XZ)  e  2 


(16) 


(17) 
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■so  19  A(I9-3) 

t  =  P  +  B(^-  +  X^)  e  1  (18) 

Substitution  of  the  above  equations  and  the  appropriate  Christoffel 
symbols  in  Eq.  (14)  gives  the  equation  of  motion  in  the  form 

*  IP  *  B(X  +  eA(1^3)I  *  |  0w  |||  <»> 


The  incompressibility  condition  leads  to: 

2  2  2  2 
*  -  Rj  =  r  -  rj 


and 


2 

tfR 

It? 


dR 


dR, 


d2R, 


i\1  ur\4  i  ur\|  rv«  u  rv1 

-5^  ‘ar*  ‘ar'  +ir  a?1 


(20) 


(21) 


With  P}»P2  denoting  the  pressure  on  the  inside  and  outside  wall, 

respectively,  of  the  blood  vessel,  the  use  of  the  boundary  conditions, 

t11  =  -pj(t)  at  R  =  Rj  and  x11  *  -p2(t)  at  R  *  R2’  substituting  Eq.  (21) 

Into  Eq.  (19)  and  integrating  yields 

d2R,  RP  dR.  2  R?  ,  R?2 

P^t)  -  p2(t)  =  pwRj  1n  jq  +  (3TJ  pw  nn  1^  +  1  (^T  *  1)1 


fX2  j  2  A(  A2  ♦  X  -  2) 

-  B  J  TT  •  *  « 


(22) 


2  2 

It  must  be  recognized  that  the  relationship  -  1  +  A  +  1/X  has  been 
used  to  obtain  Eq.  (22). 

The  following  dimensionless  quantities  are  introduced  into  Eq.  (22): 


B 


P*  -  — ,  R  *  «  ^  ,  R2*  *  ^  ,  t*  -  ,B*  *  ,o 

p  U*  1  Ro  Ko  Ko  p  u  r'o 

0  0 


*  =  —  (23) 

w  p. 
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Then  the  equation  of  motion  in  the  radial  direction  becomes 


(Pj*  -  P2*>  =  pJ  v-—?  <!q*> 


2  R  *  R  *2 

.  rt*,dR\  rin  jL  +  I  _i_  i\i 

+  PWxdt*'  D  *  +  9  R  *2  ” 


V  2  V2 


.  g*  fX2  lid  eA(x2  +  17x2  -  2)iX 

I  A  >3 


If  the  "stars”  are  dropped  for  convenience,  Eq.  (24)  can  be  written  as: 

d2R.  R?  dR.  2  Ro  ,  Rp2 

Pj(t)  -  p2(t)  «  pw«i  -r  ln  sf  -  Pw  <dtr>  nn  ^ +  ^ ' 1)1 

-  B  f  i#  eAU2  +  1/>2  '  \ 


The  initial  conditions  are: 

At  time  t  =  tQ,  Rj  =  Rq,  dRj/dt  =  u,  radial  velocity  of  fluid. 

In  the  above  derivation,  it  must  be  noted  that  only  the  radial 
displacements  of  the  blood  vessels  are  considered  significant  since  the 
axial  displacements  are  prevented  by  tethering  of  the  vessels  to  the 
surrounding  tissue. 

C.  Equations  of  Left  Ventricular  Contraction 

In  view  of  the  large  volume  changes  of  the  left  ventricle  between 
the  systole  and  the  diastole,  the  theory  of  large  elastic  deformation  is 
used  to  analyze  the  pumping  action  of  the  heart.  In  first  approximation, 
the  components  of  the  pulmonary  circulation  and  the  left  ventricle  are 
lumped  together  and  the  system  Is  treated  as  a  highly  deformable  sphere 
undergoing  radial  deformation  of  the  left  ventricle  by  using  a  strain 
energy  function  of  the  type  represented  in  Eq.  (15). 
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The  undeformed  and  deformed  configurations  of  a  spherical  chamber 
are  shown  in  Fig.  4.  Let  a  point  located  by  r,  8,  <f>  in  the  undeformed 
sphere  be  displaced  to  a  new  location  R,  0,  <|>  in  the  deformed  sphere. 

Eqs.  (12-15)  in  the  previous  section,  being  general,  are  valid  for  the 
deformation  of  the  sphere  also.  However,  the  stresses  are  expressed 
in  the  form 

=  efl<'2  -  3) 

T22.A  t  BeA(12-3,(i/x2.l4)  (», 

R4"  R2 


22 

33  _  x 

T - ~T 

sin  0 

4  2 

where  X  «  R/r  and  I2  =  X  +  2/X  .  The  equation  of  motion  in  the  radial 
direction  is 


3  (D  +  2B  A(I2-3),  +  28  A(I2-3)  ,1  .4.  9% 
^(p+^e  t  )+-Re  *  A  '  paF 


(27) 


The  incompressibility  condition  leads  to 

32R  d2Rl 


(28) 


Substituting  Eq. (28)into(27)  and  noting  the  boundary  conditions 
x11  =  -PjU) 

t22  =  -p2(t) 


(29) 


and  integrating  with  respect  to  R,  the  equation  of  motion  can  be  put 


Pj(t)  -  P2(t) 


=  .  jX2  “  e  a<a4  *  2/*2  -  3>  (X3  *  1) 


dX 


dR, 


+  o"!2  ^  <KT  '  +  2,)R1  (3T>  ^ 


<  ,J_  .  1  , 

~i  '~dt}  (rr  tt 

R1  K2 


4-) 

K2 


(30) 


Using  the  dimensionless  variables  described  in  Eq.  (23)  the  nondimensional 
form  of  the  equation  of  motion  becomes: 

Pj*  -  P2*  ■ 

fA 


-2B*  l  2  (1  +  £)  e  +  2'X  -  3> 


dX 


d2R 


?  u  *i  i  i 

+  (Rj*7p*  ^*3-  (r~t  -  flnr) 


dRl*  ,  1  1 


+  2Rj*p*  -gp-  (j^t  -  j^r) 


P*RV  A!,2,  1 

'  dt* '  p  *4  *4' 

K1  r2 


(31) 


Deleting  the  "stars"  for  convenience,  the  equation  of  motion  can  now  be 
put  in  the  form: 

Pj(t)  -  p2(t)  * 

?  fX2  n  .  1  x  _  A(X4  +  2/X2  -  3) 

-2B  J  +  KJ  >  e  < 


'dX 
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The  initial  conditions  are: 

At  time  t  =  tQ,  R^  =  RQ,  dRj/dt  =  u^,  a  time  function  which  depends  upon 
the  venous  return.  It  must  be  noted  that  the  transmural  pressure  across 
the  myocardium  in  Eq.  (32)  is  not  the  same  as  the  pressure  difference  in 
Eq.  (25). 

For  a  complete  solution  Eqs.  8,  9,  10,  25  and  32  must  be  simultaneously 
solved  with  the  appropriate  initial  and  boundary  conditions  in  conjunction 
with  a  reasonable  baroreceptor  control  mechanism. 

D.  Baroreceptor  Reflex  Control 

The  baroreceptor  control  of  the  systemic  arterial  pressure  is 

accomplished  by  a  closed  loop  regulator  which  continuously  monitors  the 

systemic  pressure  through  baroreceptors  located  in  the  carotid  sinus  and 
39 

in  the  aortic  arch  .  A  typical  steady-state  relationship  between  the 
input  pressure  (feedback)  and  the  output  pressure  (regulated  systemic 
arterial  pressure),  as  described  by  Taylor  ,  is  shown  in  Fig.  5.  Since 
it  has  been  observed  in  several  cardiovascular  system  experiments  under 
acceleration  stress  that  the  response  of  the  baroreceptor  reflex  mechanism 
begins  In  G-8  seconds  after  the  pressure  change,  it  Is  reasonable  to  use 
the  steady-state  curve  of  Fig.  5  for  model  response  under  high  sustained 
acceleration.  For  short  duration.  Impact  type  accelerations  this  curve 
would  be  unsuitable. 
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Fig.  4.  Unde formed  and  deformed  representation 
of  the  left  ventricle 


Fig.  5  Input-output  pressure  relationship  In 
baroreceptor  control 
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E.  Effect  of  Acceleration  on  Microcirculation 

The  blood  vessels  of  microcirculation  are  extraordinarily  small, 
and  their  typical  dimensions  are  of  the  order  of  microns.  Under  normal 
circumstances,  the  velocity  of  the  blood  in  the  microcirculation  is  1  mm/sec 
and  the  Reynolds  number  is  of  the  order  0(10  ),  which  is  sufficiently 
small  so  that  the  Stokes  flow  approximations  are  applicable.  Neglecting 
the  inertial  effects  and  assuming  that  the  stream  lines  are  nearly  parallel, 
the  dimensionless  equation  of  fluid  motion  in  the  axial  (2)  direction  becomes 


9w  _ 
3l 


IE 

9z 


1 


a2w  +  1 
37+  r 


9w 

9r 


2 

3w, 

3pr 


+  g(t) 


which  can  be  rearranged  to  read 


(33) 


Re !? =  - Re  <!f> +  <r! +  ?  If  +  It’  *  9<t)  l3i> 

r  z 

In  the  earth's  natural  gravitational  field,  the  dimensionless  g,  as  given 
in  Eq.  (4),  is  of  the  order  0(10~2),  and  with  the  effect  of  Re  0(10“3)  in 
the  last  term  Re  g(t)  in  Eq.  (34)  becomes  physiologically  insignificant, 
being  of  the  order  0(10"5).  We  estimate  that  the  effect  of  acceleration 
on  microcirculation  per  se  can  be  safely  neglected  up  to  100  g.  However, 
the  pressure  of  the  blood  pooled  in  the  arteries  and  veins  can  affect  the 
flow  rate  in  the  small  vessels.  For  this  reason  it  is  necessary  to 
determine  a  relationship  between  the  pressure  gradient  and  the  flow  rate 
In  the  small  blood  vessels. 

41 

For  the  flow  of  a  Newtonian  fluid  in  a  uniform  tube  Szymanski  showed 

2 

that  the  flow  would  be  fully  developed  if  vt/D  >  1,  where  t  *  time, 

v  =  kinematic  viscosity,  and  D  *  tube  diameter.  An  extension  of  this 

2 

criterion  to  microcirculation  yields  va/D  >  1  for  flow  to  be  quasi-steady, 


where  At  is  the  smallest  characteristic  time  of  the  unsteadiness  in  flow 


42 

According  to  Burton  ,  At =0.1  sec;  using  v  =  0.04  Stokes,  one  finds 

that  the  diameter  D  must  be  greater  than  600y  (microns)  for  any  significant 

effect  of  unsteadiness.  Since,  in  microcirculation  the  diameters  of  blood 

vessels  are  much  less  than  600y,  changes  in  flow  due  to  unsteadiness 

43 

become  entirely  negligible.  On  this  basis  Benis  argued  that  the  effect 
of  unsteadiness  on  non-Newtonian  flow  could  also  be  neglected.  Thus,  the 
use  of  steady-flow  equations  can  be  justified  for  microcirculation. 

For  steady  capillary  flow,  the  flow  rate  through  a  circular  tube 
can  be  expressed  by 


Q  =  2rr 


rw  dr 


(35) 


where  Q  =  flowrate,  R  3  tube  radius,  and  w  *  blood  velocity.  Integration 
by  parts  of  the  right  hand  side  yields 


Q  =  tt  |  d(r2w)  -  tt  J  r2(^)  dr 


(36) 


The  first  integral  on  the  right  hand  side  of  Eq.  (36)  is  zero.  In  the 
second  integral  the  domain  of  integration  can  be  divided  into  two  regions: 
a  cone  of  unsheared  fluid  extending  to  radius  Ry,  and  the  annular  region 
bounded  by  the  unsheared  fluid  and  the  tube  wall.  Then, 
rR..  o  j..  rR 


Q  =  -  tt  J  y  r2(|p)  dr  -  it  |  r2(~)  dr 


(37) 


The  first  term  on  the  right  hand  side  in  Eq.  (37)  which  represents  the 
core  Integral  is  zero.  By  changing  the  variables  from  r  to  t  in  the 
second  term  as  suggested  by  Merrill  et.  al.**,  Eq.  (37)  can  be  written  as 

T. 


8tt 


(ap/l) 


7 


f  w  2  •  j 

j  T  Y  dx 

'T. 


(38) 


where  L  •  length  of  the  capillary 
AP  =  pressure  drop 
t  *  shear  stress 
Y  =  shear  rate. 

The  shear  stress  and  the  shear  rate  are  related  by  an  empirical  equation 


1/2  1/2  .  1/2  *1/2 
r  '  =  t  +  u  '  v 


1/2 

in  which  is  the  yield  shear  stress,  and  y  is  a  constant  which 
represents  the  slope  of  the  Casson  plot  relating  the  viscometric  parameters 


of  blood.  In  Poiseuille's  flow  y  becomes  the  blood  viscosity.  Substitution 


of  Eq.  (39)  into  Eq.  (38)  and  Integration  yields 


■  ieW)  4*  \1/Z  R?/2  ‘tP/L)1/?  »  Tr* 

®  7/7  y  21(AP/L)3u 


it  t  R 


which  is  valid  under  the  assumption  that  the  flow  is  steady,  laminar  and 
incompressible,  and  blood  is  homogeneous.  In  the  above  equation,  and  u 
are  known  constants;  then  plots  of  Q  vs.  AP/L  for  capillaries  of  various 
radii  can  be  easily  constructed.  An  example  of  this  relationship  is  shown 
in  Fig.  6. 

RESULTS  AND  DISCUSSION 

The  system  of  equations  (8),  (9),  (10),  (15),  (25)  and  (32) 
constitute  the  basic  hydrodynamic  model  in  this  study.  These  equations 
are  nonlinear  and  coupled  and  therefore  a  closed-form  solution  cannot  be 
found.  They  were  solved  numerically  on  a  digital  computer  by  a  finite 
difference  scheme.  The  numerical  data  for  arterial  segments  were  taken 
from  Westerhof*®  and  are  presented  here  in  Table  I.  Although  the  arterial 
segments  were  presented  as  a  series  of  elastic  tubes  of  constant  radii 
which  changed  stepwise  dlstally  from  the  aorta,  a  tube  of  slight  taper  was 
considered  In  the  solution  to  avoid  the  difficulties  presented  by  discon¬ 
tinuities  in  radii. 


Table  I.  Numerical  Data  on  Arterial  Segments 


Name  of  Artery 

Length 

AX, 

(cm) 

Internal 

radius 

r 

(cm) 

Wall 

Thickness 

h 

(cm) 

Vol  ume 

Q  =  ttr2  AH 
(cm3) 

Aorta  ascendens 

2.0 

1.47 

0.164 

13.571 

Aorta  ascendens 

2.0 

1.44 

0.161 

13.028 

Arcus  aorta 

2.0 

1.12 

0.132 

7.881 

Arcus  aorta 

3.9 

1.07 

0.127 

14.027 

Aorta  thoracal is 

5.2 

0.999 

0.120 

16.303 

Aorta  thoracal is 

5.2 

0.675 

0.090 

7.443 

Aorta  thoracal  is 

5.2 

0.645 

0.087 

6.796 

Aorta  abdominal is 

5.3 

0.610 

0.084 

6.196 

Aorta  abdominal  is 

5.3 

0.580 

0.082 

5.601 

Aorta  abdominal  is 

5.3 

0.548 

0.078 

5.000 

A.iliaca  communis 

5.8 

0.368 

0.063 

2.467 

A.iliaca  externa 

5.8 

0.290 

0.055 

1.532 

A.iliaca  externa 

2.5 

0.290 

0.055 

0.660 

A. profundus 

6.3 

0.255 

0.052 

1.287 

A. profundus  femoris 

6.3 

0.186 

0.046 

0.685 

A. femoral is 

6.1 

0.270 

0.053 

1.397 

A. femoral  is 

6.1 

0.259 

0.052 

1.285 

A. femoral  is 

6.1 

0.249 

0.051 

1.188 

A. femoral  is 

6.1 

0.238 

0.050 

1.085 

A. femoral  is 

7.1 

0.225 

0.049 

1.129 

A.  pop  li  tea 

6.3 

0.213 

0.048 

0.898 

A. pop li tea 

6.3 

0.202 

0.047 

0.807 

A.poplitea 

6.3 

0.190 

0.046 

0.705 

A. tibialis  posterior 

6.7 

0.247 

0.051 

1.284 

A. tibialis  posterior 

6.7 

0.219 

0.049 

1.009 

A. tibialis  posterior 

6.7 

0.192 

0.046 

0.776 

A. tibialis  posterior 

6.7 

0.165 

0.044 

0.573 

A. tibialis  posterior 

5.3 

0.141 

0.041 

0.331 

A. tibialis  anterior 

7.5 

0.130 

0.039 

0.398 

A. tibialis  anterior 

7.5 

0.030 

0.039 

0.398 

A. tibialis  anterior 

7.5 

0.130 

0.039 

0.398 

A. tibialis  anterior 

7.5 

0.130 

0.039 

0.398 

A. tibialis  anterior 

4.3 

0.130 

0.039 

0.228 

A.anonyma 

3.4 

0.620 

0.086 

4.106 

A.subclavia 

3.4 

0.423 

0.067 

1.911 

A.subclavia 

6.8 

0.403 

0.066 

3.969 

A. axillaris 

6.1 

0.364 

0.062 

2.539 

A. axillaris 

5.6 

0.314 

0.057 

1.734 

A.brachialis 

6.3 

0.282 

0.055 

1.574 

A. brachial  is 

6.3 

0.266 

0.053 

1.400 

A.brachialis 

6.3 

0.250 

0.052 

1.237 

A.brachialis 

4.6 

0.236 

0.050 

0.804 

A.ulnaris 

6.7 

0.215 

0.049 

0.972 

A.ulnaris 

6.7 

0.203 

0.047 

0.367 

A.ulnaris 

6.7 

0.192 

0.046 

0.776 

A.ulnaris 

3.7 

0.183 

0.045 

0.389 

A.radlalis 

7.1 

0.174 

0.044 

0.675 

Table  I.  (Cont.) 


Name  of  Artery 

Length 

At 

(cm) 

Internal 

radius 

r 

(cm) 

Wall 

Thickness 

h 

(cm) 

Vol ume 

Q  =  nr2  At 
(cm3) 

A. radial is 

7.1 

0.162 

0.043 

0.585 

A.radialis 

7.1 

0.150 

0.042 

0.502 

A.radialis 

2.2 

0.142 

0.041 

0.139 

A.interossea  volaris 

7.9 

0.091 

0.028 

0.205 

A.coelica 

1.0 

0.390 

0.064 

0.478 

A.gastrica  sin. 

7.1 

0.180 

0.045 

0.723 

A.lienalis 

6.3 

0.275 

0.054 

1.497 

A.hepatica 

6.6 

0.220 

0.049 

1.003 

A. renal is 

3.2 

0.260 

0.052 

0.679 

A.mesenterica  sup. 

5.9 

0.435 

0.069 

3.507 

A.mesenterica  inf. 

5.0 

0.160 

0.043 

0.402 

A.carotis  com.  sin. 

5.9 

0.370 

0.063 

2.537 

A.carotis  com.  sin. 

5.9 

0.370 

0.063 

2.537 

A.carotis  com.  sin. 

5.9 

0.370 

0.063 

2.537 

A.carotis  com.  sin. 

3.1 

0.370 

0.063 

1.333 

A. car.  int.  sin. 

5.9 

0.177 

0.045 

0.581 

A. car.  int.  sin. 

5.9 

0.129 

0.039 

0.308 

A. cerebri  anterior  sin. 

5.9 

0.083 

0.026 

0.128 

A. car.  ext.  sin. 

5.9 

0.177 

0.045 

0.531 

A. car.  ext.  sin. 

5.9 

0.129 

0.039 

0.308 

A. car.  ext.  sin. 

5.9 

0.083 

0.026 

0.178 

A. car  com.  dextra 

5.9 

0.370 

0.063 

Zii? 

A. car.  com.  dextra 

5.9 

0.370 

0.063 

2.537 

A. car.  com.  dextra 

5.9 

0.370 

0.063 

2.537 

A. car.  ext.  dextra 

5.9 

0.177 

0.045 

0.580 

A. car.  ext.  dextra 

5.9 

0.129 

0.039 

0.308 

A. car.  ext.  dextra 

5.9 

0.083 

0.026 

0.127 

A. car.  int.  dextra 

5.9 

0.177 

0.045 

0.581 

A. car.  int.  dextra 

5.9 

0.129 

0.039 

0.308 

A. cerebri  anterior  dex. 

5.9 

0.083 

0.026 

0.127 

A. vertebral is 

7.1 

0.188 

0.046 

0.788 

A. vertebral is 

7.7 

0.183 

0.045 

0.810 

The  following  constants  were  used  In  the  solution  for  arterial  deformations: 


Rq  =  1.47  cm 
U  =11.9  cm/s 

3 

pw  =  p0  *  9r/cm 
v  =  0.038  Stoke 
A  =  0.8 

A  2 

B  =  11.35  x  10  dynes/cm 
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The  flow  chart  of  the  solution  technique  for  the  comprehensive 
model  is  presented  in  Fig.  7.  Assuming  the  pressure  variation  in  the 
heart  chamber  as  a  cosine  function  varying  between  80  and  120  ninHg,  the 
entrance  velocity  into  the  aorta  was  calculated  from  the  solution  of 
Eq.(32).  Coupling  the  linearized  Navier-Stokes  equations  and  the 
equations  of  motion  for  elastic  tube  deformation  a  solution  was  obtained 
for  pressure  and  radius  at  various  locations  downstream  from  the  aorta. 

The  volume  change  can  be  calculated  as  a  function  of  time  from  the 
computed  radii. 

Two  acceleration  profiles  were  used  in  the  solution  of  the  model 

proposed  here.  First,  an  impact  type  acceleration  which  was  experi- 

46 

mentally  determined  by  Hanson  and  shown  in  Fig.  8  was  used  and  the 
aortic  pressure  was  calculated.  This  pressure  is  shown  in  Fig.  9  in 
comparison  with  the  pressure  measured  in  the  thoracic  aorta  of  a  beagle 
dog  in  response  to  the  same  acceleration  profile.  An  exact  quantitative 
comparison  between  the  two  pressure  profiles  is  not  possible  because  they 
were  calculated  for  two  different  species  having  different  physical 
dimensions  and  material  properties.  However,  the  qualitative  comparison 
is  quite  satisfactory.  Second,  an  acceleration  profile  measured  during 
the  flight  of  an  F-14  aircraft  was  used  to  compute  both  radius  and  pressure 
in  a  segment  of  the  aorta.  The  acceleration  profile  and  the  corresponding 
radius  and  pressure  are  shown  in  Figures  10,  11  and  12,  respectively. 

The  radius  and  pressure  results  of  Figs.  11  and  12  were  obtained  for  the 
50-sec  time  domain  Indicated  In  Fig.  10.  There  are  no  pressure  and  radius 
data  available  for  the  F-14  acceleration  profile  for  comparison.  Such 
data  Is  difficult  to  obtain  nonlnvaslvely  and  there  are  understandably 
numerous  restrictions  against  invasive  measurements. 
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baroreceptor  reflex  control 

Flow  chart  for  solution  technique 


Pressure  mmHg 


Fig.  8  Experimental  deceleration  profile  measured  by  Hanson4.6 


Fig.  9  Pressure  in  the  aortic  arch  for  Hanson's 
acceleration  profile 
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Fig.  10  6-profile  for  an  F-14  flight 


Time,  sec 

Fig.  11  Radius  of  aorta  for  the  G-proflle  of  F-14 


30 


Fig.  12  Aortic  pressure  for  G-proflle  of  F-14 


It  should  be  apparent  from  Eqs.(25)  and  (32)  that  the  blood 
flow  response  will  depend  upon  the  values  assigned  to  the  material 
constants  A  and  B.  To  develop  the  total  cardiovascular  system  model 
the  values  of  A  and  B  for  all  segments  of  arteries  and  veins  in  the 
body  must  be  known.  Such  information  is  not  available  at  the  present 
time.  Also,  the  numerical  data  on  venous  segments  are  not  available 
for  the  complete  venous  system.  However,  in  the  present  study  it  has 
been  demonstrated  that  it  is  possible  to  construct  an  acceptable 
mathematical  model  of  the  cardiovascular  system  under  acceleration  stress 
by  considering  the  mechanical  parameters  such  as  the  physical  dimensions, 
fluid  properties  and  material  properties  of  the  blood  vessels. 

It  is  well  known  that  finite-difference  numerical  procedures 
are  time  consuming  and  difficult  to  formulate  for  complex  geometries 
of  realistic  shapes  of  the  left  ventricle  and  curves  and  branches  in 
arteries  and  veins.  Further  study  is  recommended  to  formulate  the  model 
by  finite  element  techniques  by  which  complex  geometries  can  be  easily 
handled  and  computations  can  be  more  economically  performed. 
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